Super Guassian Distribution of Laser Produced Plasmas 



Kaifeng Chen 



(N 

o 

(N 

> 
O 

(N 
Oh 



in 

O 

Oh 



> 
O 
(N 
On 

in 



(N 



i 



Contents 

1 Introduction U 

2 Physical Model U 

3 Simulated Results H 

4 Conclusion 

5 References [U 



2 



1 Introduction 



It is well known that the non-Maxwellian electron velocity distributions, resulting from the 
nonclassical drive and transport under laser-fusion relevant conditions, may significantly af- 
fect the dispersion and damping of the different waves present in the plasma. This in turn 
may alter the thresholds and gains of the parametric instabilities, involving the plasma wave 
and the ion acoustic wave, respectively. 

A systematic characterization of the electron velocity distribution, under the wide range of 
experimentally relevant conditions, is therefore essential in this context. Electron distributions 
have been computed and commented in many previous publications. But, cither because they 
considered idealistic profiles, or considered only a small number of simulation cases, these 
studies did not attempt to characterize the non-Maxwellian distribution by proposing, for 
instance, a functional form as an appropriate fit. 

A thorough study of the non-Maxwellian electron velocity distribution functions was per- 
formed by Mattee el.al., neglecting spatial transport. In this studythe competition between 
inverse Bremsstrahlung, which pushes the distribution towards a super-Guassian of the 
form/ ~ exp[— (v/u) 5 ],&nd the thermalization through self-collisions was considered.lt was 
shown, that the distribution is well represented, at all energies, by so so-called Dum-Langdon- 
Mattc(DLM)-type function Jdlm ~ exp[—(v/u) n ], where the parameter n takes intermedi- 
ate values between n = 2 (a Maxwellian),in the absence of drive, and n — 5, when inverse 
Bremsstrahlung dominates. As a practical result, a relation n = n(a) was established be- 
tween DLM exponent n and the relevant parameter a = Z(v os /vth) 2 , measuring the relative 
importance of inverse Bresstrahlung compared to self-collisions. Here,Z is the degree of 
ionization, v os is the quiver velocity and v t u the thermal velocity. 



2 Physical Model 

The transport of electrons is described by the corresponding Fokker-Planck equation: 

where E stands for the self-consistent electrostatic field ensuring quasineutrality,C ee stands 
for the Landau electron-electron collision operator, and C e i for the Lorentz electron-ion 
collision operator. As for the simulations, the ions are kept fixed. 

For solving the Fokker-Planck equation, a high ionic charge Z is assumed, such that the 
electron-ion mean free path \ e i is small compared to any characteristic scale length of the 
system. So the distribution function can be a sum of an expansion series. 

oo 

f(x,v,t) = J2fi(x,v,t)P l ( f i) (2) 

1=0 

where fx — v x /v = cos9 is the direction cosine. In the diffusion approximation. 

A factor a = Z(v os /v t h) 2 is usually introduced to describe the relative strength of these two 
processes. When a is small, the thermalization process is dominant so that electron have a 
distribution close to a Maxwellian. If the heating process is dominant, i.e. a ids significantly 
larger than 1, the electron distribution should deviate from a Maxwellian and approach to 
a super-Guassian distribution, 

f DL M(v,t) = A n ^exp[-B n { — )»] (3) 

v th V th 
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where the coefficients A n and B n , ensuring the definition of density (AT = 4tt J °° v 2 dvf) and 
temperature(3 NT /m — Att v 4 dvf), are given in terms of the gamma function T by 



A,. 



ir(5/n) 



B n = 



47i-r(3/n) L3T(3/n) 
1 T(5/r 



3/2 



</2 



-3T(3/n) 

Here m is a number depending on the factor a, 

m(a) -2 + 3/(1 + 1.66/a - 724 ) 



and v m is given by 



t k B TT(3/m) 



(4) 
(5) 

(6) 
(7) 



m e r(5/m) 

In this work, the data is generated by Doc. Bin Zhao's Fokker-Planck code in Epper- 
lein's scheme. And these data contain information about the distribution function of 
(x, v, t), temperature of each point and density of each point. All these data are fitted by 
using the formula proposed above. 



3 Simulated Results 




Figure 1: x=75 

Figure 1 shows the densitytemperature and distribution function at the 75th lattice. As to 
the distribution function, the x axis has been transformed to u 2 , which represents the energy 
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distribution, and the y axis denote the logarithm of f. Thus for Maxwellian distribution 
function like / = Ciexp(—C2V 2 ), it will be a straight line and the slope will be 2(represented 
in the green line). The blue line is plotted by using DLM fitting function and the red line 
denotes the results calculated by Fokker-Planck equation. Some parameters are listed above 
each sub-picture like m,the ratio of n e to n^and the temperature(-£T). 

In this picture, the ratio is 1.55, which means that the density is much beyond the critical 
density. And for this region of this high ratio, m — 2 and the Fokker-Planck calculated data 
are little deviated from the DLM theory. 

Figure 2 presents the lattice of x = 76, the ratio becomes nearly 1.27 and the three lines: 
pure Maxwell distribution, super-Guassian distribution function and the distribution function 
from Fokker-Planck equation. In this position, these three lines are almost the same because 
m = 2. 09, which is very close to 2. 




Figure 2: x=76 

Figure 3 shows some slight changes when the position becomes x — 77. In this case, the 
distribution calculated by Fokker-Planck equation is less than the super-Guassian fitting 
function at high velocity region. More over, these two functions are always higher than the 
Maxwellian distribution function. This phenomenon agrees well with the non-local transport 
theory that in high velocity region of high density heat flux is actually lower than the classical 
one, thus the percentage of high velocity particles will become larger. 

Under this condition,™ = 2.11 and the ratio of n e /n c is nearly 1.00. The Figure 4 shows 
the results when n e < n c , which means that this position along x axis is the point away 
from the over density region. In this region, the Maxwellian distribution function becomes 
the largest one and the DLM function and the Fokker-Planck calculated result are nearly 

1 The critical density. 
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Figure 3: x=77 




Figure 4: x=78 



G 



the same. More to mention, from about v = 3v t h the DLM and Fokkcr-Planck begin to 
differ from each other. This phenomenon proves the validity of Brunner's fitting function in 
reference 1. The value of m now is 2.15. 

4 Conclusion 

In summary, in the case of spacial transport and inverse Bremsstrahlung, some simulations 
by using Fokker-Planck equation are shown. Once the system has reached a quasisteady 
state, the electron velocity distributions within the laser beam are well represented by a 
DLM- type function for the bulk. By considering the ratio of n e and n c ,we can better 
assume the distribution functions. In the high density region, Maxwellian distribution can be 
used like a classical model. However, when the density goes below the critical density, DLM 
function should be used to simulate the non-local transport. Some evidence has shown that 
when v is close to 2>v t h, there will be some deviation of the DLM- type function and Fokkcr- 
Planck calculated data. Future work will contain other cases and continuing to explore 
better fitting functions in the case of inverse Bremsstrahlung. 
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